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The RCSJ model of resistively and capacitively shunted Josephson junctions is used to describe 
superconducting point contacts over a wide range of resistances up to the metallic - tunneling 
transition. Their small dynamic capacitance of order C — 0.1 fF due to the point-contact geometry 
results in a huge plasma frequency. The critical current is then strongly suppressed and the contact 
resistance becomes finite because of quantum-mechanical zero-point fluctuations of the Josephson 
plasma and the rather large escape rate out of the zero-voltage state due to quantum tunneling. We 
test the predictions of the RCSJ model on the classical superconductors lead, indium, aluminum, 
and cadmium. 



I. OVERVIEW 

Point contacts between two superconductors have attracted interest for a number of reasons: They can be used to 
study the behaviour of a macroscopic quantum objectu as well as quantum transport when the junction consists of few 
atoms onlyErtl. They may also help to extract information on the order parameter of compounds like the heayy-fermion 
superconductors, although those experiments are not as easy to interpret as has previously been assumedu. n n 

Metallic point contacts of classical superconductors have first been investigated systematically by Muller et al. uu. 
They observed that, on reducing the contact area, the residual contact resistance Rq at zero bias becomes finite and 
the critical current I c smaller than the theoretical value I® at a temperature T -C T c even when the normal-state 
resistance i?N = R{T > T c ) is far below the quantum limit Rk/2 = h/2e 2 — 12.9 kil. This was attributed to the 
presence of electrical radio-frequency noise and additional damping of the junctions at high frequencies due to the 
small impedance (that is the large capacitance) of the current and voltage leadsQtl. The Andreev-reflection excess 
current, on the other hand, was found not to depend on the lateral size of the junctions. 

We show here that the behaviour of those Josephson junctions can better be described by the small capacitance 
C typical for the three-dimensional point-contact geometry. Taking into account the capacitance of the junctions as 
an adjustable parameter, we derive relations between i?N and the residual resistance Rq as well as the product R.^I C , 
and estimate both the plasma frequency cu p and the Josephson coupling energy Ejk from the current-voltage I(U) 
characteristic at I — > 0. Our experiments on the classical superconductors lead, indium, aluminum, and cadmium 
fit the predictions fairly well. Using the horizon model, the capacitance C can be traced back to the properties of 
vacuum tunneling junctions in the normal state. 

II. THEORY 

The system of Cooper pairs in a bulk sample fornas one wavefunction \t = |^o| e-xp (i</?). The gradient of the phase 
ip drives the current. According to B. D. Josephsoncl, by coupling two superconductors weakly to each other using a 
thin insulating layer between them, the phase difference ip := <f2 — </?i across the junction results in a supercurrent 

J = 7° sin p (1) 

as long as the voltage drop U is much smaller than the superconducting energy gap 2A(T)/e. According to Ambe- 
gaokar and BaratoflO the critical current of a Josephson tunnel junction = (7rA/2ei?N) tanh (A/2fcsT). 

The critical current through short metallic junctions between isotropic BGS-type superconductors was derived by 
Kulik and Omel'yanchuk to be (KOI) J° = (1.32 nA/2eR N ) tanh (A/2fc B T)El in the dirty limit (electron mean-free 
path I much shorter than the superconducting coherence length £), and (K02) 1° = (7rA/ei?N) tanh (A/2fceT)E£l in 
the clean limit I 3> £. We will see below that quantum fluctuations of the Josephson plasma can further reduce the 
critical current. The superscript denotes then the above intrinsic value. 
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The current-phase relationship of a metallic Josephson junction in tlie clean limit deviates from the simple sinusoidal 
of a tunnel junction. It depends on the transmission coefficient D astJ 
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with S = yj\ — D sin 2 (<p/2). The normal-state contact resistance in the ballistic limit 
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Here fcp is the Fermi wave number and a the contact radius. By varying the transmission coefficient, the junction can 
be tuned from pure metallic (D = 1) to pure tunnel (D — 0) behaviour. 

Note the effects we are going to investigate do not involve an 'energy cost' due to lateral confinement that could 
suppress the superconducting features in the contact region: although the junctions are rather small - with a typical 
metallic fcp ~ 15 nm -1 a i?N = 100 SI junction has a contact radius a of about 2 nm - the Heisenberg uncertainty 
principle yields a spread with respect to momentum, and not one with respect to energy. 
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FIG. 1. RCSJ model. The junction itself is represented by x. The current is applied to and the voltage detected at the 
terminals denoted by o. 

A. The RCSJ model and the quasi-particle resistance 

When the applied voltage U ^ 0, the phase does not change because <p = dip/dt = 2eU/h — 0, and the current 
through the junction is constant!! If the voltage is finite, the phase increases at a rate 2eUlh ~ 3 TO 15 (1/s) -U{VX and 
the current has an oscillatory component at extremely high frequencies. D. E. McCumberO and W. C. Stewartta took 
into account that besides the supercurrent across the Josephson tunnel junction there exist also a normal quasi-particle 
current, described by a parallel resistance R — R qp , as well as a displacement current due to a (stray) capacitance 
C (Fig. |l|). This capacitance is important for the properties of the junction because it represents a short circuit for 
the high-frequency part of the supercurrent. The total current of such a resistively and capacitively shunted junction 
(RCSJ) is I = 1° sin <p + U/R qp + CdU/dt or 
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with the plasma frequency uj p = ^/2eI®/hC and the characteristic frequency lo c = 2eI®R qp /h. 

The quasi-particle resistance R = R qp is a priori unknown. Generally, it will depend both on voltage, on current, 
or on phase. Above T c , without Josephson effect, the quasi-particle resistance equals just the normal-state contact 
resistance _Rn, originating purely from the quasi-particle current through the junction (we assume there is no external 
shunt resistance). For a tunnel junction, this R qp is proportional to the square modulus of the tunnel matrix element. 
For a ballistic junction, like the ones we are using, R qp involves an infinite resummation of multiple quasi-particle 
transmissions through the contact. Then, energy is not dissipated at the junction itself: electrons equilibriate far away 
from the junction in the bulk metal. By including non-linear effects like quasi-particle quasi-particle interactions or 
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electron-phonon coupling, the contact resistance, and thus R qp , depends on temperature and the applied bias voltage 
U. 

Below T c , the situation changes completely as a Josephson supercurrent can flow through the junction. However, 
well below T c there is no quasi-particle density of states within the energy gap 2A on both sides of the junction. 
Therefore, R qp — > oo for an ideal junction at low bias voltage \U\ <C A/e. At a Josephson tunnel junction, R qp can 
be directly observed by suppressing the supercurrent using a small magnetic field less than the critical field. Tunnel 
experiments have shown that the damping resistance can be much mailer than the intrinsic quasi-particle resistance 
i? qp , see for example Ref.EJ. This is probably due to parasitic damping of the high-frequency electromagnetic field 
generated at the junction. 

The quasi-particle resistance of a metallic Josephson junction (or point contact) can not be measured directly, as 
far as we know. But it seems plausible, to apply the same arguing as for the tunnel junctions: For a clean BCS-type 
superconductor far below T Cl the quasi-particle density of states has vanished in the vicinity of the chemical potential, 
and therefore the intrinsic R qp — > oo. Recently, Levy Yeyati et allB derived the dissipative part of the current through 
a superconducting quantum point contact, that is a ballistic junction with one conducting channel. Assuming the 
same relation to hold also in the general case of a ballistic Josephson junction, the invers quasi-particle resistance 
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r = %/t being the inelastic relaxation rate, and r the quasi-particle lifetime. Most importantly, the dissipative part 
of the current vanishes at <p — > or at / <C I c , that is R qp ^> i?N ■ 

In both type of experiments, Josephson tunnel junctions as well as metallic Josephson point contacts, there can 
be damping due to the external circuitry. The amount of this parasitic damping can not be predicted theoretically, 
because it depends on the specific experimental setup. Thus the damping resistance R = R qp of the RCSJ model has 
to be treated as a free parameter, like the capacitance, that has to be extracted from the properties of the junctions. 
The normal-state contact resistance may represent some lower bound of the intrinsic quasi-particle resistance, and 
i?N niay thus serve as a first guess for the damping resistance. But otherwise there is no reason to set R qp = R^. 
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FIG. 2. Washboard potential E(I, ip) of Eq. H for different values of the applied current. 



Eq. U describes the motion of a particle with 'mass' (h/2e) 2 C exposed to a viscous damping force (h/2eR qp )(p in 
the potential 



E(I, tp) = - —lip - E m ■ cos<^ 
2e 



(6) 



Eje = ^J<?/2e is called Josephson coupling energy. The potential Eq. || has the shape of a washboard, its wells have a 
depth of 2_Eje at I = (Fig. ^). The number of discrete energy levels in these wells amounts to aboutEllnt(2i?jE/?iWp). 
The lowest one is the zero-point energy eo = hw p /2. Minima exist as long as I < I®, and the static solution of Eq. [| 
is I /I® = sin if. If the time constant R qp C is much smaller than the period of the plasma oscillations, the system is 
overdamped. The only stable state is then at such a potential minimum and the I(U) - characteristic is non-hysteretic. 
An underdamped junction can cither be at a minimum, the supercurrent flowing without a voltage drop. Or it can 
run down the washboard, and the current is accompanied by a voltage. The I(U) - characteristic is then hysteretic. 
At I > 7° there are no more minima, and the system is definitely in the resistive state with ip =/= 0. 

The RCSJ model has been applied to weakly damped planar tunnel junctions (with an insulatiag oxide layer instead 
of a vacuum gap) of classical superconductors like jaiobium and tin: Fulton and EUinklebergerlij observed thermally 
activated escape from the minima, Voss and WebbliU and later on Washburn et al&3 reported macroscopic quantum 
tunneling through the wells of the washboard potential. Its discrete energy levels have been observed by Martinis et 
al.iB. Obviously, the RCSJ model is the minimum system to study the Josephson effect at solitary point contacts. 
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FIG. 3. Josephson coupling energy -Bje, zero-point energy hu} p /2 at two different capacitances C, and thermal energy kBT 
at 50 mK vs. Rn/T c of junctions between BCS-type superconductors in the clean limit. 

Let us illustrate the situation by setting in numbers. For junctions with BCS-type superconductors in the clean 
limit the coupling energy amounts to 



E JE w 1 eV 



Tc(K) 



at a zero-point energy of the Josephson plasma 



12 meV , 



To(K) 



fl N (n)c(£P) 



(7) 



(8) 



Planar tunnel junctions typically have capacitances of order 1 pF. For them hw p /2 can be orders of magnitude smaller 
than Eje- To induce a detectable escape from the minima of the potential wells, these junctions have to be driven 
near the critical current. Metallic point contacts, on the other hand, have considerably smaller capacitances. They 
are believed to be of order 1 fFEl. And although the coupling energy can be much larger than the thermal energy 
fceT (about 4/ieV at our lowest temperature of ~ 50 mK), the zero-point energy can be quite high, see Fig. ||. This 
has direct consequences with respect to the critical current and the residual resistance, because in such a system 
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the particle can escape easily from the minima even at small currents, accompanied by 2tt phase slips. In addition, 
because of the quality factor 



Q = uj p R qp C « 0.04fl qp (fl) ^ c( ^^ fF) (9) 

those metallic junctions have to be expected to be in the cross-over region between weak and strong damping, if the 
quasi-particle resistance i? qp was of the same size as the normal-state resistance i?N- 



B. Critical current suppressed by quantum fluctuations 



The quantum-mechanical treatment of the particle in the washboard potential describes the phase as an operator 
conjugated to the charge q — 2e of the Cooper pairs, [tp, q] = 2ie. This allows to include quantum fluctuations and 
phase diffusion. The macroscopic voltage is then given by the time derivative of the expectation value of the phase 
U = (H/2e)d (<p(t)) jdt. It was shown, for example, that Josephson junctions at T = can undergo a phase transition 
from the superconducting to the resistive state, depending on the strength of dissipation or dampingc3. _ . 

A current 7 = xl® tilts the washboard potential (Fig. |2|) and reduces both the minimum height of the wcllsEJ 



2E(x) = Eje [x{2 arcsin x — n) + 2 cosarcsinx] 



(10) 



and the plasma frequency 
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7 C can closely approximate J° at small i?N when fkj p <C -Eje- The 7(f ^characteristic shows then a sudden rise at 
7 as 7,?. At larger resistances, Eje decreases faster than fkJ p , and the ratio Hoj(x)/E(x) diverges as x — ► 1, that is 
at I — > I®. Assuming the supercurrent being suppressed by the quantum- mechanical fluctuations of the Josephson 
plasma, the actual critical current I c of the RCSJ model less than 7° is reached at 



2E(x c ) » huo(x c )/2 



(12) 



when the lowest energy level coincides with the minimum height of the potential well. The phase difference across the 
junction is then not well-defined any more. Eq. p"2| makes the reduced critical current x c = 7 c /7° an implicit function 
of the intrinsic 7^ and of the capacitance via 
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Because of the quantum fluctuations the 7(f) - characteristic is washed out and also the critical current becomes less 
well defined. 

Eq. O relates the normal-state resistance of a junction to its reduced critical current x c as 



i? N (x c ) «i? K (i? N 7 c °) 



C 



lu p E(x c 
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(14) 



depending only on the capacitance C and the material-dependent parameter (R^I®)- The supercurrent vanishes 
completely above 



R N {x c = 0) w Rk{RnIc) 



C 



(15) 



Although we have neglected damping, Eq. [15] resembles the s = y '^R^CI® / 'en =1.5 and r] = 1 
fixpoint for the localization - derealization transition at T — > 0, the so-called Schmid transitior 
that our simple physical argument captures the essence of a renomalisation group calculation. 
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FIG. 4. Energy of the stationary states in units of 2_Ejb vs. reduced mass C* = Eje/2Ec for tunnel junctions (D = 0) and 
metallic junctions (D — 1) at low (T/T c — 0.01) and high (T/T c = 0.1) temperatures (from top to bottom). The analytical 
zero-point energy eo = 2_Eje \J -Ec/2-E.je as well as the first and second excited level (solid lines) fit well the numerical data for 
the tunnel junctions (Ec = e 2 /2C). 

A metallic Josephson junction with finite transmission coefficient slightly changes the situation with respect to a 
tunnel junction. Neglecting dissipation, the Schrodinger equation 

can be solved numerically. Fig. || shows the Eigen energy e = q of the stationary states i = 0,1,2,... fit quite well 
the analytical data of the washboard potential Eq. |^ at small energies. For metallic junctions with the current-phase 
relationship of Eq. ^, the potential minima broaden and the levels shift towards smaller energies. Since these are small 
corrections, for ease of use, and because we do not know the transmission coefficient of the individual junctions, our 
discussion will be based on the washboard potential. 

For i?N^c ~ 0.1 meV and C » 1 £F the critical resistance of Eq. |l5| is of order 5 kS~l, which is still in the metallic 
regime. It shows that these phenomena due to the capacitance of the junctions have to be taken into account seriously. 



C. Phase diffusion: Finite contact resistance and suppression of the critical current 

At sufficiently high temperatures thermal activation induces escape from the potential minima by 2tt phase slips, 
that is the phase diffuses. In this regime, the finite contact resistance due to a small capacitance has already been 
discussed by Ambegaokar and HalperinEJ. 

At T — * escape fiam-the potential wells is mainly due to quantum tunneling. The tunneling rate has the form, 
see for example Refs.BEj, 

r QT = ^exp(-£) . (17) 
In WKB approximation and neglecting damping, the coefficients 



G 



A = A Q 




and 



D = - / y/C[E(tp)-eo]dtp 



(18) 



(19) 



The iptegral runs over that phase space in which E(ip) > €q. The parameter Aq = V60 for the cubic potential of 
Refs. E3'E3. (We did no succeed in deriving Aq, but we hope that it does not depend strongly on the shape of the 
potential). We have estimated the coefficients A and B both for tunnel and metallic junctions (Fig. ||). Neglecting 
the exact shape of the washboard potential, both A and B may be overestimated by about 20 per cent. In the 
tunneling rate these corrections partly compensate. Our numerical results strongly deviate from the analytical form 
when the mass becomes small. The zero-point energy is then quite high, the effective potential the particle has to 
tunnel through becomes small, and the WKB approximation breaks down. 
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FIG. 5. a) A/Aq and b) B vs. reduced mass C* = Em/2Ec for tunnel junctions (D = 0) and metallic junctions (D = 1) at 
low (T/Tc = 0.01) and high (T/T c = 0.1) temperatures. The analytical A/A = 2e y / B/2n/h and B = ^/WEje/Ec (solid 
lines) fit well the numerical data of the tunnel junctions at large masses. 



Caldeira and LeggettEl and Larkin et aZ@ derived the rate 



QT 



7 -exp 



14.4E JE 

TlLJ-o 



0.87 



(20) 



with 7 w -\/1207r(14.4£ , jE/?2.Wp) for tunneling out of a cubic potential. Eq. |20| holds exactly at weak damping, when 
the quality factor Q> 1. The tunneling rate is considerably reduced at strong damping Q <^ 1. In what follows, and 
for ease of calculation, we assume the junctions are only weakly damped by including the term linear in l/Q. 

Note that for the sinusoidal current-phase relation and the washboard potential (that is for / <C I c ) the prefactor 
in the exponent is not 14.4 as for the cubic potential (that is for I rs J c ) but 8y/2 w 11.3. Since this is also a small 
correction, like for the A and the B coefficients mentioned above, we will discuss our experimental results using the 
above tunneling rate Eq. GO. This also has the advantage of including possible corrections due to damping. 
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FIG. 6. Residual resistance Rq vs. _Eje/^ p at weak damping Q ^> 1, Q = 5, and Q = 1. 

The asymmetry of the washboard potential due to a current flow yields a net phase slip rate ip, thatpis the phase 
drifts. Tilting the potential enhances the tunneling rate exponentially. Following Tinkham's discussionEJ we get for 
the RCSJ model and Eq. 20 at weak damping 



ip(x) ~ 47iTqt sinh 



7.2%xE JE 

htOr, 



0.87 



(21) 



Here we have used E(x) sa -Eje(1 — irx/2) and u>(x) w u> p , a good approximation as long as x — I/I® < 0.5. The 
phase slips result in a voltage drop U — htp/2e. Consequently the differential resistance 



dU (I) I dl « R cosh ( —7 



0.87 



(22) 



with the residual resistance Rq = dU /dl(l = 0) 



Ro » I.&jRk 



0.87 



exp 



TlUJ-n 



0.87 



(23) 



At Q 1 one can directly read off the Ejft/hui p -ra,tio from Ro in Fig. ||. At strong damping, the phase becomes more 
localized and Rq being strongly reduced. The second-order approximation of Eq. |22] 
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(24) 



allows to derive the plasma frequency uj p and, through i? , the coupling energy _Eje from the I(U) - characteristic 
at / <C ij. In this approximation neglecting damping means to underestimate the plasma frequency by a factor of 
[1 + 0.87/Q]. The derivation of the Josephson coupling energy, on the other hand, is barely affected by damping 
because the correction factors almost cancel each other. 

To ensure tunneling from the lowest oscillator level at hu p /2, the tilt of the washboard potential hi /2e per period 
must be less than the energy difference of about huj p between the two lowest levels. At J > Iz = eujp/ir quantum 
tunneling populates the second level, enhancing the total tunneling rate. This may be called Zener tunneling of the 



8 



phase. And because the phase diffuses faster than described by Eq. ^D, the differential resistance should increase 
much stronger stronger than predicted by Eq. |2^. 

The supercurrent is suppressed as soon as the phase diffuses or drifts faster than w p or at 



1. 



This relates the critical current 
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i? K < Ro sinh f ^-1 
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(26) 



with the residual resistance and the plasma frequency. Here we have assumed that the Q-factor does neither depend 
on / nor on U. At Q > 1 the critical current is larger than Iz as long as Rq < 1.8Rk/ sinh(3.6) « 2.6 kf2. Eq. |26| 
represents an upper limit for I c because of neglecting Zener tunneling. Note that this equation is valid only when the 
actual critical current has been reduced already, that is at x c < 0.5, while Eq. ^ represents a maximum for all x. 

An estimate like in the previous section for the critical current yields for a i?N = 1 junction a coupling energy of 
Ejb ~ 0.2 meV and a plasma frequency of uj p « 0.57 (ps _1 ). Neglecting damping, the residual resistance due to phase 
diffusion Eq. ^3] amounts then to Rq ss 0.92 kO. Because of the exponential dependence in Eq. a i?N = 100 SI 
junction has a tiny Rq ss 58 only. Thus the residual resistance rises steeply within a very narrow resistance range. 



D. Bloch-wave oscillations 



To add a charge q to the capacitance C of the junction requires the Coulomb charging energy Eq — q 2 /2C. When 
Ec > Eje, both phase and charge do no longer behave as classical variables, but like operators. Because of the 
periodic washboard potential (Eq. ^|) the Josephson junction catt then be described using Bloch wavesafia. Their 
wave number is the quasi charge k = q/2e. For a review see RefB3. 

Bloch-wave oscillations represent the periodic transfer of discrete Cooper pairs across the junction, which is 
recharged by the external current source at a radial frequency of ttI /e. These wave packages travelling in (p space of 
the tilted washboard potential have a finite 'momentum' < dip/dt >. This implies a voltage drop, and the contact 
resistance becomes finite. TJaeocy predicts a region at finite bias current in which Bloch waves contribute negatively 
to the differential resistancd±3'E3, but the details of the spectra are difficult to derive unless the junctions are in the 
Eq ^S> Eje limit. Since Bloch waves are suppressed by fluctuations, the quasi-particle resistance, that is the resistance 
R qp of the RCSJ model, must be larger than about -Rk/4 ~ 6.4 kO. It is therefore unlikely to observe them at 
metallic junctions if the quasi-particle resistance had the magnitude of i?N- Intraband transitions can be induced 
when the external current is strong enough. This Zener tunneling of the quasi-charge requires / > ew p /7r, like for 
Zener tunneling of the phase, see the previous section. _ 

There are only few experimental reports on Bloch-wave oscillations at Josephson junctions, e. g. Ref.cZl. In those 
experiments Bloch oscillations have been identified by a systematic coincidence between structures of the I(U) - 
characteristic at certain applied DC currents and the frequency of an external microwave excitation. 



E. Andreev reflection 



Additional information on the properties of the junction comes from Andreev reflection at finite voltages. If 
the potential difference is \U\ > A/e, quasi-particles from one side of the junction can enter the other one by 
Andreev reflection, and return a quasi-hole. At high bias voltages this additional hole o r e xc ess current amounts to 
I cx = 8A/3ei?N and I cx = (tt 2 — 4) A/4ei?N in the clean and the dirty limit, respectivelyc§Ea. 

If \U\ < A/e the quasi-hole returns at an energy still inside the gap, and it is reflected again as a quasi-particle. 
Since the charge carriers gain an energy e\U\ per crossing the junction, M = 2A/e|Z7| (M = 1,2,3,...) successive 
reflections are required to overcome the gap. At \U\ — 2A/eM the conductance increases stepwise. This multiple 
Andreev reflection becomes more pronounced at the presence of normal quasi-particle reflection at the interface, that 
is at D < 1. Normal quasi-particle reflection strongly reduces the excess current, because both the quasi-particle and 
the quasi-hole have to cross the interface. 

The Andrccv-rcflection excess current can also be reduced due to a finite quasi-particle lifetime r. For normal- 
superconducting junctions, lifetime effects are taken into account by the modified BTK theoryEJ, applying the Dynes' 
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model. At superconducting junctions, lifetime effects are expected to suppress structures due to multiple Andreev 
reflection, see for example Refs.Hu. 

As soon as phase diffusion becomes efficient and Rq finite, it can be difficult to distinguish clearly between the 
Josephson supercurrent and the Andreev-reflection hole current. At large voltages \U\ 3> A/e and beyond the critical 
current, however, the Andreev-reflection process is not affected by the capacity of the RCSJ model. Therefore ±he 
literature value of J ex has to be expected. This has been verified previously on normal-superconducting junctionsEj. 
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FIG. 7. Our break-junction device. The sample is broken by turning the micrometer screw using two flexible cotton threads. 
This screw is also used for coarsely adjusting the contact. Fine adjustment is achieved by applying a voltage at the piezo tube. 




FIG. 8. The sample on the bending beam. This gold-plated beam is cut from copper alloy. Only one of the current and one 
of the voltage leads are shown. 
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III. EXPERIMENT 



We prepared metallic point contacts between two bulk pieces of the BCS-type superconductors lead (T c = 7.2 Kh 
indium (3.4 K), aluminum (1.15 K), and cadmium (0.54 K) using mechanical-controllable break junctions. Fig. [7] 
shows our setup and Fig. || the sample on the bending beam. The samples, ~ 1 mm diameter wires, were broken at 
a predefined notch in the ultra-high vacuum region of the cold refrigerator. This avoids oxidation of the interfaces 
and ensures the junctions to consist of pure metal. The lateral contact size and, thus, the normal-state resistance 
could be adjusted in situ by a micrometer screw, driven mechanically by a pulley-and-rope and a piezo tube. The 
current-voltage characteristic and the differential resistance were recorded in the standard four-terminal mode with 
current biasing. 

Twisted pairs of wires are used for the current and voltage leads. They are properly shielded. Simple LRC filters 
at the mixing chamber protect the sample from low-frequency (~MHz) noise. We have no copper-powder filter as 
described in Refsflia. 

We start the experiments with the contact resistance set to about 1 — 3 mfl at room temperature. This corresponds 
to a contact radius of about 10 /*m. The residual resistance ratio of the material in the contact region was then 
observed if possible while cooling down. We obtained as lower bound of the electronic mean free path / > 70 nm 
(Pb), 310 nm (In), 545 nm (Al), and 75 nm (Cd). The junctions-investigated here have i?N > 1 f2. They are therefore 
ballistic and have, according to Sharvin's resistance formulatEl, that is Eq. g at D = 1, a contact area of about 
665nm 2 /i? N (f7). The bulk samples had residual resistance ratios of about 7570 (Pb), 4760 (In), 1270 (Al), and 1430 
(Cd). This corresponds to electronic mean free paths of about 18 fim (Pb), 29 /im (In), 19 fim (Al), and 12 fxm (Cd). 
Unless otherwise specified, the following experimental data have been obtained at low temperatures of T w 50 mK. 
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FIG. 9. I vs. U and dU/dl vs. U of a) a low-resistance and b) a high-resistance In junction at T = 50mK(thick solid lines). 
Arrows define the various parameters used. The vertical bars mark the superconducting energy gap 2A/e and its subharmonics. 
The dashed lines are the normal-state I(U) characteristics. The thin solid lines represent the additional current due to the 
Josephson effect and the Andreev-reflection excess current. 



IV. RESULTS AND DISCUSSION 



The I(U) - characteristics as well as the differential resistance dU/dl could be recorded over the full superconducting 
anomaly without excessive heating for junctions with i?N > 1 f2. Few junctions with i?N ~ 1 — 10f2 were hysteretic, 
either due to heating or, according to the RCSJ model, due to being weakly damped. The self- magnetic field at the 
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junctions could also contribute to hysteresis, but this has to be expected at much smaller contact resistances i?N < 1H- 
For those hysteretic junctions the critical current was extracted from that branch of the I(U) - characteristic with 
increasing current. 
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FIG. 10. Typical spectra dU/dl vs. U of a) Pb and In and b) Al and Cd junctions at T — 50 mK. 



Fig. g shows two typical I(U) - characteristics and dU/dl vs. U spectra as example and for defining the various 
parameters. A number of selected spectra are shown in Fig. [l^. They look quite similar for all four superconductors. 
There are distinct structures due to superconductivity at small i?N, including an unresolvably small R . These 
structures diminish at larger resistances, but their positions on the voltage axis do barely change. Finally, in the 
resistance range around 10 kil we observe, as expected, the transition to the tunneling regime, with zero-bias maxima 
instead of minima. The contact resistance i?N up to which well pronounced superconducting anomalies could be seen 
was the higher the larger the superconducting gap was. 

The gap structure at U = 2A/e and its subharmonics due to multiple Andreev reflection agreed well with the BCS 
value 2Abcs = 3.52 fceTc = 1.05 meV, 352 jUeV, and 165 /zeV of In, Al, and Cd, respectively (Fig. pd|) . Strong-coupling 
Pb had the width enhanced by about 25 per cent with respect to 2Abcs = 2.20 meV (Fig. |ll|). Why the Pb junctions 
did not have the right energy gap at small resistances we do not know yet. Heating effects or the suppression of 
superconductivity by the self-magnetic field can be excluded because the other superconductors behave regularly (at 
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least for junctions with i?N > 10 ft). Some of the Al junctions had lower T c down to 0.85 K, while some of the Cd 
junctions had higher T c up to 0.65 K. This is probably due to disorder and stress in the contact region and resulted 
in smaller or larger gaps, respectively. 
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FIG. 11. Experimental superconducting gap 2A at T = 50 mK vs. 7?n, normalized to the BCS value 2Abcs = 3.52 /cbTc. 
Identical symbols are used in Figs, [n] - ^ for the same pieces of sample. 
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FIG. 12. Experimental 7? N / e x vs. i? N at T = 50 mK, normalized to the cleandimit (K02) value 8A/3e = 3.67, 1.40, 0.47, 
and 0.22 mV for Pb, In, Al, and Cd, respectively. Solid lines indicate the K02 (= 1.0) and KOI (= 0.55) value. 



A. Andreev-reflection excess current 



At small i?N ~ 1 ft, the excess current I ex was near the clean (K02)limit for Pb (using the experimental 2A = 
2.75 meV) and In, and near the dirty (KOI) limit for Al and Cd (Fig. |lj). This trend could be consistent with the 
superconducting coherence length £ oc 1/A being shorter for Pb and In than for Al and Cd. However, according to 
the above bulk resistivity data all four superconductors should have been in the clean limit. 
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At larger i?N and towards the metallic-tunneling transition at -Rk/2, the excess current of the Pb and In junctions 
decreased but remained finite, settling around the dirty limit. This reduction may have several reasons: First, small 
contacts could be stronger distorted than large ones. Second, theJunctions could have a finite quasi-particle density 
of states due to lifetime effects. Using the modified BTK theorvo, a reduction of the excess current by a factor two 
requires a lifetime parameter not larger than r = Ti/t ss A/5. It agrees roughly with r ss 50 /ieV, estimated using the 
bulk residual resistivity of the samples. This estimate would be supported by the clearly resolved multiple Andreev 
reflection anomalies of Pb and In, while Al and Cd show more washed out and smeared structures. Third, normal 
quasi-particle reflection. Since multiple Andreev reflection is observed, the transmission coefficient D has to be at 
least somewhat less than 1. An excess current of (on the average) not less than half the maximum possible one, 
makes the lower bound of the transmission coefficient D > 0.5. On the other hand, all four superconductors show the 
transition from vacuum tunneling to metallic conduction at around i?x/2, indicating D»l. 




R N (£2) 



FIG. 13. Experimental R^I C vs. _Rn at T = 50mK, normalized to the clean-limit (K02) value 7rA/e = 4.32, 1.65, 0.55, and 
0.26 mV for Pb, In, Al, and Cd, respectively. Solid lines indicate the K02 (= 1.0) and the KOI (= 0.66) result as well as that 
due to the zero-point fluctuations (Eq. [uj) at C = 0.05 fF and the abscissa fixed at 0.8, the best fit for aluminum. 



B. Critical current 



At small i?N « 10, the critical current I c (Fig. [T^) was near the clean (K02) limit for Pb and In, and near the 
dirty (KOI) limit for Al and Cd, like the excess current. But R^I C decreased continously at larger contact resistances 
and vanished well below -Rk/2. This agrees with previous observations by otheraJu. 

Like for the excess current, we would have expected the critical current or the product R^I C to be reduced, but to 
the dirty limit only. Assuming a minimum D sa 0.5 for normal reflection from the Andreev-reflection excess current 
and applying the current-phase relationship of Eq. ^ the critical current can be reduced by almost a factor two 
near the Ambegaokar-Baratoff value 7rA/2ei?N- Thus the reduction of R^I C to zero requires an explanation that 
does neither affect the superconducting gap nor the excess current. The RCSJ model could solve this problem: the 



high-resistance part of R^I C can be described by Eq. 14 and a small capacitance of C w 0.05 fF. The fit is nearly 



perfect for Al, but deviates for the other superconductors. Especially the nearly straight line of Pb in the semi-log 
plot astonishes. These deviations from the theoretical curve at constant C require additional processes, for example 
a variation of the intrinsic R^I® a t the- junction, a varying amount dissipation, or a varying capacitance. 
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FIG. 14. Experimental -Rn^c vs. Ro at T = 50mK, normalized to the clean-limit (K02) value 7rA/e. Indicated are the K02 
(= 1.0) and the KOI (= 0.66) result as well as that due to phase diffusion Eqs. ^ and |4l| at a lead capacitance of k — 3.3pF/m 
and a quality factor Q S> 1 (solid lines) and Q — 1 (dashed lines). 
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FIG. 15. Residual resistance Ro vs. Rn at T = 50 mK. Straight lines indicate the detection limit, dashed lines are Ro of 
Eq. ^ at weak damping Q ^> 1, using the clean-limit 7° and C = 0.05 fF. Solid lines take into account Rn/2 in parallel. 

Fig. |l4| shows an alternative plot of R^I C vs. residual resistance Ro- Again a systematic decrease is observed when 
Rq increases. But this time, all four superconductors behave similarly. The fit through the data points takes into 
account the dynamic capacitance of the junctions as discussed later. 
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C. Residual contact resistance 



At low i?N the residual contact resistance Rq = dU/dI(I = 0) was unresolvably small. Here the experimental 
resolution results from the ~ 2 nV detection limit at a maximum current of I c . Towards larger i?N the residual 
resistance emerged out of the detection limit of ~ 2nV/J c , first rising steeply, arj|d then approaching ~ Rn/2 (Fig. fTo] ). 
Again, this behaviour agrees qualitatively with previous observations by othersETEl. 

Raising the temperature, for example to T — 2K for indium, did barely alter the results. This indicates that the 
junctions are not thermally activated but are in the regime of quantum tunneling. Only near T c , when the intrinsic 
(i?N^° ) decreases strongly, we also found the transition to a finite Rq at smaller (Fig. 

The data in Fig. |l5|, especially the steep rise of Rq (Rn), can well be described by Eq. |2^ using a constant C w 0.05 fF 
as for the critical current in the previous section, the clean-limit I®, and weak damping. Deviations at large i?N have 
to be expected: first, the above analysis does not apply when the experimental RnI c = 0. Second, in the dissipative 
state the total resistance of the junction can not exceed i?N- It must be less than Rn/2 if one includes Andreev 
reflection, that has not been considered by Eq. |2^. Combining phenomenologically these two different regimes, that 
is Rq of Eq. ^3] and Rn/2, by substituting 1/Rq — > 1/Rq + 2/R^, fits the data of all four superconductors. 
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FIG. 16. a) Experimental RnIc vs. _Rn of In at the indicated temperatures. The data are normalized to the clean-limit 
(K02) value 7rA/e at T — > 0. Solid lines result from the zero-point fluctuations Eq. [H] at C — 0.05 fF. b) Residual resistance 
i?o vs. i?N- Temperatures and symbols as in a). The straight line represents the detection limit for the T — 50 mK data, solid 
lines are Rq of Eq. G3| using the clean-limit I^iT) and C = 0.05 fF, and taking into account the asymptotic -Rn/2. 



D. Plasma frequency and Josephson coupling energy 

According to Eqs. ^ - 12^ the plasma frequency as well as the coupling energy of each contact with a sufficiently large 
residual resistance can be derived from the differential resistance at small currents. We find indeed a linear dependence 
between dU/dl and I 2 (Fig. ^). Its slope divided by Rq yields directly the plasma frequency u> p , assuming weak 
damping Q 1. This in turn is used to extract -Eje from Rq. The contacts that can be analyzed this way have 
residual resistances of ~ 0.1 fl to ~ 1 kfi, and according to Fig. ^ their Ej^/huj p varies from about 1.2 to 0.5. 



16 




1/1° 

FIG. 17. Typical dU/dl — Ro vs. 1/1° at T = 50 mK. Here 7° is the clean-limit value derived from 7?n. The dashed lines 
have a slope of 2 in the double-log plot. 
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FIG. 18. Josephson coupling energy Eje vs. Rn derived from dU/dl at 7 — > and at T — 50 mK, assuming weak damping 
Q ^> 1. Solid lines are Eje ~ A7?k/47?n in the clean limit. 

In contrast to the strong reduction of the experimental I c in Fig. |l3|, the coupling energy 7?je at I — * in Fig. [IJ is 
slightly larger than expected from the clean-limit 7°. The plasma frequency u> p corresponds to a capacitance of order 
C = 0.05 fF only at large contact resistances (Fig. [l9|). Deviations at small 7?n can not be explaind by a reduced 
escape rate due to stronger damping, because including damping would lead to an even larger lo p and not to a smaller 
one. Thus we are forced to conclude that the capacitance C is not a constant. The effective capacitance of the 
junctions 
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C^-P? (27) 

becomes C ss 0.1 fF/7°(/xA), depending on the intrinsic critical current (Fig. ^). C does not depend on the geometrical 
cross sectional area of the junctions which is inversely proportional to R^. We will show below that it is the above 
relationship with I® that has indeed to be expected. 
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FIG. 20. Capacitance of the junctions C derived from Eq. ^ and the data of Figs, [[i] and |l^ vs. J?n- Solid lines represent 
C = 0.1 fF/J°(^A) and the clean-limit J c °. 

For the junctions with finite Rq, the -EjE/?iWp-ratio is always close to one, that is these junctions are very susceptible 
to quantum fluctuations and phase diffusion. And it requires only slight variations of this ratio to produce the steep 
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rise of Ro(R^) because of the large prefactor 14.4 in the exponential function of Eq. 

If one assumes i? qp = J?n, the damping factors Q = 8tt R qp EjE / RkJv^p ~ ^qp(k^) may justify our data analysis 
using the Q 3> 1 limit for the tunneling rate Fqt and neglecting damping for Pb junctions at Rjsi > 1 kf2 but not 
for Cd, despite both superconductors show similar results. The question arises whether it is correct to describe the 
damping term of these ballistic junctions by their normal-state resistance i?N- Our experiments seem to indicate that 
the true quasi-particle resistance of these ballistic Josephson junctions at I -C I~ is much larger than Rjq. Such an 
enhanced i? qp has to be expected according to Eq. |], it finds its natural explanation in the gap of the quasi-particle 
density of states, that is there are no quasi-particles at all (the clearly resolved multiple Andreev-reflection signal shows 
that lifetime effects are negligible). Damping can then only be due to losses of the high-frequency electromagnetic 
field around the junction in the sample and in the metallic part of the setup. 

Simultaneously, the absence of hysteretic /([/^-characteristics that would be typical for weakly damped junctions 
can be explained by strong phase diffusion or drift at J c that hinders trapping of the phase. One may also speculate 
whether the Q factor near I c is smaller than that at J — ► 0, due to the frequency dependence of the electromagnetic 
field generated by the alternating Josephson supercurrent. 




I(|iA) 

FIG. 21. dU/dl vs. I at T = 50 mK of a Pb and an In junction. The dashed lines are fits using Eq. ^| The lower two 
diagrams show the experimental data normalized to the fit. Arrows mark the current I* below / c at which the measured curves 
starts to deviate from the fit. 



E. Zener tunneling of the phase 

-Eje and w p have been derived from the spectra at rather small currents. The intermediate range below the 
experimental I c usually shows a steep rise of dU/dl. Some of the junctions, especially those with Pb and In, even 
have a distinct kink at I* w 0.2/° (Fig. ^l|). Such a steep rise has to be expected, and it is straightforwardly explained 
as Zener tunneling of the phase at /* — Iz, when the lowest energy level of one of the potential wells exceeds the 
next higher level of the neighbouring well. This is the more important since Eje ~ hujp of these junctions implies 
that the washboard potential contains not more than two or three discrete levels. Fig. E2 shows that the agreement 
between I* and Iz is quite good for large uj p (and small Rq). At small frequencies (and large Rq), the current I* 
falls below etOp/n. This has to be expected, too, because the level spacing of the anharmonic washboard potential 
decreases when Ej-E/hu> p decreases. 
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FIG. 22. I* vs. plasma frequency ui p . The solid line is 7* = Iz- The squares mark the anomalies attributed to Zener 
tunneling of the phase, circles denote the position of the minima of the spectra. 

Nevertheless, Zener tunneling of the phase implies that the junctions are only weakly damped. Otherwise, the 
experimental Rq would yield E^/tiuj-p < 0.5, and the washboard potential would have not more than one discrete 
level. As Zener tunneling features could be observed at junctions with Rq up to about 260, 175, and 15 O for Pb, In, 
and Al, respectively, the quality factor Q > 5. This supports our interpretation in the previous section. 

Zener tunneling of the phase is preferably resolved at Pb or In junctions and not at Al or Cd junctions. The simplest 
explanation is that the voltage drop at I = Iz 

sinh(3.6) „ 

O.D7T 

of about 0.26 fiV ■ u; p (l/ps)i2o(fi) has to be considerably smaller than the superconducting gap 2A/e. This condition 
is the easier fulfilled the larger the energy gap is, that is for Pb and not for Cd. 
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F. Possible excitation of Bloch waves 



One source of uncertainty in deriving both Eje and w p are additional structures - zero-bias maxima - superposed 



on dU{I)/dI of some of the junctions (Fig. 23). They were observed quite often for Al and Cd and less often for 



Pb and In. Conventional explanations for such zero-bias maxima are Kondo scattering at magnetic impurities or 
Andreev reflection at a normal-superconducting interface. Both explanations seem to be feasible here, because they 
only require a small normal interfacial region or a tiny amount of Kondo impurities (on the ppm level), that may have 
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been introduced by cutting the nut. But these two phenomena should take place at nearly constant voltage, while 
the experimental data are clearly related with a current of I*/I% ~ 0.1. 

An alternative explanation is Coulomb blockade with its negative contribution to dll/dl, and one-jmay speculate 
whether these additional anomalies in Fig. ^ result from the excitation of Bloch waves, like in Ref.Ell. In favour of 
this interpretation is that in the previous paragraph we have found hints that the damping is smaller than expected 
and the quasi-particle resistance much larger than i?N- However, our junctions have Coulomb charging energies 
Ec = (hojp) 2 /8-Eje ~ -Eje/10, that means they are in the tight-binding limit. If there was a Coulomb blockade then 
it can not be expected to be well pronounced, like what is seen in Fig. [23|. The transport of Cooper pairs could also be 
enhanced when the periodic transfer of Cooper pairs is in resonance with the zero-point oscillations of the Josephson 
plasma. The position I* of the minima marks then the Zener current Iz for tunneling of the quasi-charge. As a cross- 
check, we estimate -Eje from by assuming the theoretical I® and calculate then uj p from the experimental R . 
Indeed, the position /* of the minima fits reasonably well the above estimated plasma frequency, that is ttI* /e ~ uj p , 



see Fig. 22 



A negative contribution could be superposed on the dU / dl spectra of the other junctions, even if a minimum cannot 
be resolved. The relative size of such a contribution may depend on damping, that is on the (unknown) quasi-particle 
resistance R qp . The plasma frequency and, consequently, the coupling energy would then have been overestimated. 
However, in view of the available experimental data, these possible corrections are difficult to handle. 



G. Noise 



Up to now we have neglected electrical noise. Since noise can also drive the junctions towards the resistive state, 
we have to ensure that it is small enough. A natural measure for the noise magnitude is the superconducting gap 2 A. 
Low-frequency noise up to about 2 MHz makes no problem, we can detect it directly: including amplifier noise, it is 
less than about 2/ieV (10/xeV) at a 100 SI (lOkfi) resistor and comparable to the thermal noise k^T ss 4/zeV at our 
lowest temperature of 50 mK. This is much smaller than 2A of the four superconductors, and it is consistent with 
the fine structure of the spectra due to multiple Andreev reflection in Fig. |l(]. 

The possible radio-frequency component of the external noise, however, can only be revealed indirectly through the 
properties of the junctions. For this reason we can not prove that this sort of noise is negligible. Let us assume, for 
the moment, that the finite Rq as well as the reduced I c mainly result from external noise. What is then the minimum 
noise level required to explain these results ? Weakly damped junctions in the thermally activated regime are used for 
this estimate because including damping and assuming quantum tunnelling would require more noise. The tunneling 
rate is then given by the well-known Arrhenius law 



r TA =-^exp -- — ^ (29) 
27r V k B T cfi J 

T e g is an effective temperature T c g = T + Tnoisc, that is we assume white noise with an equivalent temperature 
T no j S e, and superpose it on the thermal noise of the junctions. Using the same arguments as for quantum tunneling 
Eqs. pb]-p4| the differential resistance at I — * 

dU(I)/dI » R cosh ( * h I) (30) 



2efcBT ff 



Here the residual resistance R® = dll/dl(l = 0) 



'4fc B T cff V k B T eS 

The second-order approximation of Eq. p0| 



^K^expf-^L ) 



^-l + if-j!* ] (32) 



R dl 2 V2efc B T eff 

now allows to derive the effective temperature of the junctions from the spectra at I — > 0. By comparing the 
Independence of Eq. |32] with that of Eq. ^J, we can directly read off 

kBT e ff = ftwp/3.6 (33) 
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from the data in Fig. (the u> p in Eqs. and |33| denote different parameters). This yields lower bounds for the 
hypothetical noise temperature in the 1 - 100 K range. If such a noise level was real, we would never have seen any 
superconducting features of our samples. Thus we conclude that the external noise level is much smaller than the 
zero-point energy of the junctions, and ui p in Fig. ^ is the plasma frequency. 

The shape of the /(^-characteristics at I c responds more sensitively to external noise than the residual resistance. 
The pronounced kink at the critical current of an ideal /([/^-characteristic can be rounded, and both the hysteresis 
andpL can be reduced when the noise level exceeds even a small fraction of the coupling energy -Eje, see for example 
Ref.Ej. To observe this noise, its size must be at least comparable to the zero-point energy of the particle in the tilted 
washboard potential. Since increasing the current means to reduce the minimum height of the potential wells 2E as 
well as the zero-point energy fiuj/2, at some current below the intrinsic critical one, thermal or external noise should 
definitely overtake. However, in contrast to the analytical Eq. [H], the zero-point energy has a lower bound as described 
by Eqs. [l^ and |l3[ This is due to the fact that the flow of supercurrent needs at least one discrete energy level. Those 
junctions with finite Rq, for which we can derive the plasma frequency, would require noise temperatures in the 1 - 
1000 K range. Such a noise level is by far too large, as found above while discussing the possible effects of noise on R 
and the /(£/)-characteristic at I — ► 0. Towards smaller contact resistances, the experimental EjE/htOp ratio increases 
slightly, the critical current approaches the theoretical I® (at least for Pb and In), and the minimum w can be much 
smaller than w p . But the absolute value of cj p also rises if we extrapolate the data in Fig. |l^, so the necessary noise 
level will not go down. The hysteresis found at some of the low-resistance junctions in the i?N = 1 — 10 17 range 
could be explained by low damping, although heating effects can not be exluded. At larger contact resistances, the 
minimum required noise level becomes smaller than about T no j SC rj IK. It is then quite possible that external noise 
contributes to the spectra. We conclude that the rounding of the /(iJ)-characteristics as well as the reduction of I c 
mainly origines from the quantum-mechanical zero-point fluctuations as an internal noise source of our Josephson 
junctions and not from external noise. Of course, it is this internal noise that can also strongly reduce the hysteresis 
of the I(U) -characteristics. 

There are additional arguments that noise can indeed be neglected: 

a) Our experimental observations, agree qualitatively with previous ones by others who used a similar type of setup 
for preparing the break junctionsuu. If external noise was the limiting factor, then the noise level should be the same 
in both experiments. But it is unlikely to have an identical noise level at these various experiments and at different 
times and locations. 

b) A different series of break-junction experiments on doped semiconductors (Germanium) using our setup indicated 
a noise level of less then about 50 for junctions with i?N < 100 kfi, and corresponding to a noise temperature of 
less than about 0.5 K. Thus the total noise at the junctions consist mainly of the low-frequency component. 

c) If the high-frequency noise would dominate, it would roughly be related to the experimental zero-point energy 
huj v /2, as discussed above. We expect junctions with larger i?N to pick up more noise then those with smaller i?N, 
independent of the superconducting gap. Quite contrary, the experimental lo p in Fig. 19 decreases almost like A//?n- 
Note that for these junctions ftw p /2 > 2 A. 

d) Both the experimental critical current I c from I(U) at large currents and the I(U) characteristic at / — > 
consistently indicate a small capacitance of order 0.1 fF. 

e) The experimental Eje derived from dU (I) / dl at / — > nearly coincides with the theoretical coupling energy 
(Fig. |lj). There is no adjustable parameter. 

f ) Although rather different processes (Zener tunneling of the phase as well as the excitation of Bloch waves and 
Zcncr tunneling of the quasi-charge) may contribute to the spectra at I* ~ /z, this marker closely corresponding to 
eajp/ir supports our derivation of the plasma frequency. 

We believe these arguments as a whole strongly support the basic concepts of our interpretation. What is needed 
now is to understand the physical meaning of the capacitance of metallic Josephson junctions. 



H. The horizon model of tunnel junctions 



Future research on superconductors with even lower T c , those that have gaps with nodes of zeroes, or gaps with 
unknown symmetry like the heavy-fermion superconductors, requires to reduce the quantum fluctuations by increasing 
the capacitance. 

It is difficult to extract the capacitance of solitary junctions. This_prxiblern has been discussed extensively—for the 
case of Coulomb blockade at tunnel junctions, see for example Refs.E3~E3, applying Nazarov's horizon modelEZl. Our 
break junctions in the vacuum-tunneling regime, when the two halves of the sample are not in direct contact but 
are separated by a 0.1 — 1 nm wide vacuum gap, and in the normal state have an asymptotic offset of their I(U) 
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characteristic at ~ 50 mV in the range 0.2 - 2 mV. Attributing this offset to Coulomb blockade, we estimate an 
intrinsic capacitance C'o of order 0.1 fF due to the vacuum gap. 
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FIG. 24. a) Normalized conductance of a tunnel junction with a constant 1/7? = dl/dU (dashed line), modified by Coulomb 
blockade Eq. Bq (solid line). The relative size SR/R of the zero-bias anomaly determines the lead capacitance k, assumed to be 
n = 5pF/m. The capacitance Co, here 0.1 fF, is inversely proportional to both the width SU and SR/R. b) Experimental data 
of Al junctions in the normal state at T — 50 mK. Superconductivity has been suppressed by a magnetic field of B — 100 mT. 
7? is the differential resistance at U ~ 100 mV. The two lower spectra have been shifted for clearity. The solid lines are fits 
using Eq. ^ with k ~ 5 — 6pF/m. The intrinsic capacitance Co ~ 20 aF. 

At smaller voltages the lead capacitance enhances the total capacitance with respect to Co- The leads - which 
means here not only the current and voltage leads but also the two halves of the sample and the close surroundings of 
the junction itself - can be treated as a transmission line with a capacitance per length k. The Heisenberg uncertainty 
principle now defines the shortest time scale for the relevant processes at the junction, and thus a maximum spread 
or horizon of aboutEH ch/e\U\ (at U = 1 mV this spread amounts to about 188 /jm). Here c denotes the speed of light. 
The total capacitance becomes then 



C{U) « C + 



nch 

W\ 



(34) 



diverging at U — ► 0. To incorporate Coulomb blockade into the low-temperature spectrum of a tunnel junction with 
constant conductance 1/R = dI(U)/dU, its I(U) characteristic is displaced by e/2C(U) on the voltage axis, that is 



RI^U- 



sign(f/) 



The normalized conductance becomes 



2C(U) 

^dl ^ e 2 Kch 
dU ~ 2 (e\U\C + Kch) 2 



(35) 



(36) 



At large voltages dl/dU = 1/R is recovered, but there is a SR/R — e 2 /2nch deep mimimum at U = 0, see Fig. |24 
With this dip the lead capacitance amounts to 



2ch SR 



(37) 
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Co can then be estimated from the width SU of the anomalies. For our vacuum tunnel junctions (Fig. gfj) we get an 
average k ~ 5pF/m (for the different junctions k was found to vary between about 2 pF/m and 10 pF/m). This is a 
quite reasonable value, although being-smaller than typical literature data reported for tunnel junctions, for example 
k w 10 pf/mEa and k ps 20 - 32 pF/rr£l 
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co p (1/ps) 

FIG. 25. Capacitance C = 4e 2 -EjE/?i 2 Wp vs. plasma frequency uj p of the Josephson junctions. The solid lines represent 
C = k,c/u p with k = 3.3pF/m. 

I. The dynamic capacitance of metallic Josephson junctions 

Metallic junctions have an intrinsic Co = 0, that is only the dynamic capacitance contributes. The small experimen- 
tal capacitances in Fig. ^0| together with the lead capacitance K estimated in the previous section reveal an horizon 
of C/k ~ 1 — 100 /mi, restricted to the immediate vicinity of the junctions. The huge capacitances (~ 0.1 nF) of 
the external current and voltage leads do not matter, in contrast to the interpretation in Refs.Qfl. How to derive 
theoretically the horizon of a ballistic Josephson junction ? 

The theory for tunnel junctionsEHEa is based on the idea that the charge transfer across the junction creates an excess 
charge on the other lead. This excess charge is transported down the leads via photons with arbitrary low energy. 
The new state is orthogonal to the equilibrium ground state, yielding an infrared divergency. In the quantum regime 
of the junction, the divergency manifests itself as a zero-bias anomaly with a power law. We do not know whether 
such a theory can be applied to metallic Josephon junctions investigated here. However, we will use the physical idea 
that the infrared divergency, caused by the interaction of the system with its electromagnetic environment, is cut off 
at the lowest characteristic frequeny of the problem. 

The RCSJ model maps the many degrees of freedom of the electron system involved in charge transfer across a 
Josephson junction onto a single degree of freedom, the phase difference ip between the two bulk superconductors. 
Since phase and charge density are canonically conjugated variables, and the electromagnetic gauge field couples to 
the charge, its quantum fluctuations modify dynamically the capacitance of the junction. Moreover, the interactions 
with the electromagnetic field are strong enough to define the capacitance (or mass) of the particle in the washboard 
potential. The ground state of the quantum-mechanical model of a Josephson junction even at finite current I < I c 
is a meta-stable one. It can be described by the ground state of an harmonic oscillator with the Josephson plasma 
frequency uj p and a zero-point energy tku p /2. A constant expectation value < <p > leads to U = 0, although the charge 
fluctuations of the ground state take place on a time scale l/w p . Since this is the only time scale of the system, we 
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expect ui p to play the crucial role as cut-off frequency in denning the horizon. And this will yield a self-consistent 
picture for the interpretation of our experimental results. 

A different approach is known in quantum mechanics: a conventional particle of mass m interacting with a relativistic 
quantum field. This interaction includes a certain environment around the particle, limited by a cut-off parameter 
~ ch/mc^3. This cut-off - which is just the horizon we are looking for - depends on the internal energy mc 2 of 
the particle. It serves to avoid the creation of additional particles, excluded in the non-relativistic treatment of the 
particle itself. Applying this model to a Josephson junction lacking a real mass, the internal energy mc 2 has to be 
replaced by the excitation energy fuv p to avoid transitions between the ground state and the next higher level. 

Both models thus predict the same horizon ~ c/u> p . The total capacitance becomes then 

C re kc/uj p {C) (38) 

Note that to p is a function of C with a positive feedback: the higher the frequency the smaller is the capacitance, 
further enhancing the plasma frequency. Solving for C and lu p yields 



and 



2e/° , x 

* ^ (40) 

respectively. The experimental data in Figs. |l^ and ^o| reproduce indeed such a dependence from 1° cx l/i?N- The 
absolute values of C(u> p ) of the four investigated superconductors in Fig. 25 fit almost perfectly a k ~ 3.3pF/m. The 
good agreement with k pa 5pF/m estimated above for the tunnel junctions, supports our interpretation. 

Since frequency and capacitance are related to each other, Eqs. ^ and|4(] are valid only at I — ► 0, like Eq. |27|. When 
the applied current approaches the intrinsic I®, the capacitance should become larger because the plasma frequency 
decreases. Near J° the capacitance should diverge (this is partly compensated by the applied current introducing 
an additional time scale or frequency ~ For this reason the relation Eq. [14] can not adequately describe the 

experimental data (disregard the correct order of magnitude of C), that is quantum fluctuations only indirectly 
suppress the critical current. 

Most of the critical current data in Fig. |TJ have been reduced with respect to the theoretical 7°. Therefore it seems 
reasonable to assume C(I C ) ~ C(I — > 0). Phase diffusion at \<p\ > u> p becomes then responsible for I c being too small. 
As the Zener current Iz = euj p /ir = 41® / kcRk, the critical current at Q 3> 1 amounts to 



7cRi o^^ arcsm H^rJ (41) 

This I c depends directly on the residual resistance. Eq. ^l] fits the experimental data quite well, see Fig. [l4|. No 
adjustable parameter is required because the lead capacitance k is known from the spectra at I — > 0. Considering 
damping reduces the critical current of Eq. |l| by a factor of about [1 + 0.87/Q]. The close coincidence between 
RnIc(Rq) and Eq. ^ thus excludes any strong effect from damping. It supports our above estimate of Q — uj p R qp C = 
KcR qp > 1, and the quasi-particle resistance of the ballistic junctions being considerably larger than R^. 

There seem to be several possibilities to reduce the lead capacitance k and drive the junctions towards the delocalizcd 
Coulomb blockade limit, or to make it larger and get more localized Josephson-like behaviour. To enlarge the lead 
capacitance k by increasing the accessible volume or by bringing the junctions closer to the conducting-ground plane 
than is possible with bulk samples, one should use junctions made of thin films like that described in RefO or whiskers. 
On the other hand, the conduction electrons around the junction of a bulk sample reduce the effective volume of the 
capacitor, as long as ui p is smaller— than the plasma frequency ~ 10 4 ps _1 of the conduction electron system, as 
already proposed by K. K. Likharevc3. In an extreme limit, the Josephson junction is completely embedded in normal 
metal. Alternating currents and fluctuations could then exist only at frequencies above the plasma frequency of the 
conduction-electron system, enhancing the quantum fluctuations by driving the effective capacitance towards zero. 
Such junctions with strongly reduced critical currents may have been realized already with the so-called heavy-fermion 
superconductors. At junctions with these superconductors a considerable part of the contact region seems to be driven 
normal due to stress and disorder, and the product R^I C was found to be several orders of magnitude smaller than 
expected, see e. g. Ref.Q. To take into account a small C may be an interesting new aspect for the future study on 
those superconductors. 
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V. SUMMARY 



The quantum-mechanical treatment of the RCSJ model explains quite well the properties of superconducting point 
contacts over a wide range of resistances. The small capacitance of the junctions reproduces both the suppression of the 
critical current and the finite contact resistance by phase diffusion and drift due to quantum tunneling at the presence 
of large quantum-mechanical fluctuations of the Josephson plasma. Our interpretation does not require external 
noise and the huge capacitances of the current and voltage leads. A detailed analysis of the I(U) characteristics and 
dll(l)/dl spectra reveals a frequency-dependent capacitance of the junctions that is well described by the horizon 
model, that has first been applied to tunnel junctions. The lead capacitance is the only free parameter with which one 
can almost fully understand our Josephson junctions. Combining the properties of the metallic Josephson junctions 
with that of a completely different type of vacuum tunneling junctions in the normal state strongly supports our 
interpretation. There are several hints that the metallic Josephson junctions possess a rather large quasi-particle 
resistance i? qp > 1 kO. They are therefore only weakly damped, and anomalies which we have attributed tentatively 
to the excitation of Bloch waves, Zener tunneling of the quasi-charge, and Zener tunneling of the phase become 
plausible even though the normal-state resistance of the junctions is less than J?k/4. 
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